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Effect  of  lateral  boundary  conditions  on  electron  states  in 
tunnel-coupled  quantum  wells 

O.  E.  Raichev  and  F T.  Vasko 

Institute  of  Semiconductor  Physics,  NAS  Ukraine,  Pr.Nauki  45,  Kiev,  Ukraine 

Abstract.  The  electron  transport  in  tunnel-coupled,  double-layer  nanostructures  of  a finite 
size  is  studied  theoretically.  The  boundary  conditions  for  matrix  distribution  function  of  the 
electrons  near  edges  and  contacts  are  derived.  The  quantum  size  effects,  associated  with  the 
tunnel  coherence  length,  are  demonstrated  for  several  kinds  of  planar,  laterally-restricted  double- 
layer systems. 


Introduction 


The  fabrication  of  nanostructures  based  upon  two-dimensional  (2D)  electron  layers 
assumes  lateral  restriction  of  2D  electron  gas  as  well  as  creation  of  contacts  to  it.  Thus, 
two  types  of  boundaries,  the  edges  and  the  contacts,  are  present  in  a finite-size  2D 
system.  An  investigation  of  electrical  properties  of  such  a system  requires  a knowledge 
of  boundary  conditions  for  both  kinds  of  boundaries  and  applications  of  the  boundary 
conditions  to  description  of  the  distribution  function  of  the  electrons. 

In  this  paper  we  derive  the  boundary  conditions  for  double-layer  2D  electron  systems 
[1],  A considerable  experimental  and  theoretical  interest  to  these  systems  is  concerned 
with  manifestations  of  an  extra  degree  of  freedom  due  to  the  tunnel  coupling  of  electron 
states.  Mathematically,  it  means  that  the  Hamiltonian  of  the  system  is  represented  by 
2x2  matrix  in  the  left  (/)  and  right  (r)  layer  basis.  As  a consequence,  the  distribution 
function  of  the  electrons  is  also  a 2 x 2 matrix,  and  the  boundary  conditions  for  this  func- 
tion should  be  of  a matrix  form.  Below  we  obtain  the  boundaiy  conditions  for  columnar 
wave  functions  at  contacts  and  edges.  Then,  we  derive  the  boundaiy  conditions  for  the 
matrix  Wigner  distribution  function.  Finally,  we  apply  these  boundaiy  conditions  to 
calculations  of  the  electric  current  along  a narrow  strip  [2]  formed  from  a double-layer 
2D  system  and  across  narrow  double-layer  regions  formed  by  independently  contacted 
2D  layers  [3], 


I Wave  functions  of  finite-size  double-layer  systems 


Consider  a double-layer  system  confined  in  y direction.  The  ffee-particle  Hamiltonian 
of  the  system  is  written  as 


h1  d 2 d 2 

^W  + ^ + h- 


(1) 


Here  px  is  the  electron  momentum  along  OA,  and  h = (A/2)<rz  + Tax  is  the  potential 
energy  matrix  expressed  through  the  Pauli  matrices  &■„  A is  the  splitting  energy  in  the 
absence  of  tunnel  coupling,  and  T is  the  tunneling  matrix  element.  The  two-component 
wave  function  of  the  Hamiltonian  (1)  is  given  by 


*00  = G+[Aieik+y  + A2e~ik+y\  + g.[B^k-y  +B2e~ik-y}.  (2) 
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Here  Q±  are  the  columns  determined  by  parameters  A and  T,  and  = [2m  (e  =f 
At/2)  - pl]/h2  (A r = v/A2  + Jffy  is  the  splitting  energy  of  the  tunnel-coupled 
states).  The  columns  Q±  describe  symmetrized  G+  and  antisymmetrized  Q _ electronic 
states,  and  'F  ( y ) contains  a linear  combination  of  these  two  states. 

The  coefficients  A\,  A2,  B\,  and  B2  are  to  be  found  from  boundary  conditions.  For 
a boundary  at  y = yo  they  are  given  by  the  matrix  boundary  condition  of  third  kind: 

= 0.  (3) 

y=yo 

Here  V is  a diagonal  matrix:  V = qil)  + q,Pr,  where  I)  = (1  + <r-)/2  and  P,  = 
(1  - <j,)/2  are  the  projection  matrices.  The  complex  constants  qi  and  qr,  characterizing 
the  boundary,  depend  on  the  electron  momentum.  They  are  different  from  each  other, 
because  the  boundary  potentials  for  different  layers  are  different.  An  explicit  expression 
of  qi%r  for  a model  of  hard-wall  edge  potentials  is  presented  in  [2].  According  to  Eq.(3), 
the  total  microscopic  current  through  the  boundary  is  equal  to  (ih/2m)^'+(yQ)(V  - 
P+)'F(yo)-  Therefore,  if  the  boundary  is  an  edge  (no  current  through  it),  both  qj  and 
c/r  are  real.  If  the  boundary  is  a contact,  q\  and/or  qr  contain  imaginary  parts  proportional 
to  the  microscopic  currents  injected  in  proper  layer.  In  practice,  the  contact  can  be  made 
independently  to  / or  r layer.  If  only  / (or  /•)  layer  is  contacted,  we  have  1m  qr  = 0 (or 
lm  qj  = 0).  Different  kinds  of  structures  can  be  realized  experimentally  by  means  of 
independent  contacting  [1,4];  two  of  them  are  shown  in  Fig.  1 (see  also  [3]). 


dy 


*G0 


2 Boundary  conditions  for  Wigner  distribution  function 

The  density  matrix  of  electrons  obeys  the  quantum  kinetic  equation  dp/dt+(i/h)[H  ,p}=  0, 
where  H is  the  Hamiltonian  of  the  system.  If  one  uses  the  two-coordinate  representation 
of  the  density  matrix,  p = pPx(yuy2),  the  boundaiy  conditions  for  it  directly  follow 
from  Eq.(3),  provided  that  the  motion  of  electrons  across  the  edge  or  contact  area  is 
nearly  ballistic: 


_d_ 

dy\ 


Ppx{y\,y2)  + 'PpPs{y\,y2) 


= 0, 


J y i =>;o 


g j^Pps  (yuyi)  + pPx  (y\ , y2)V " 


= 0. 


1 y'2—yo 


(4) 


For  a finite-size  system  whose  width  is  large  in  comparison  with  the  quantum  length 
(i.e.  with  the  wavelength  of  an  electron  at  the  Fermi  level),  it  is  more  convenient 
to  introduce  Wigner  representation  of  the  density  matrix.  Then,  the  quantum  kinetic 
equation  reduces  to  a quasi-classical  kinetic  equation  for  Wigner  distribution  function 
fPxpy(y),  provided  that  this  function  changes  on  the  lengths  larger  then  the  quantum 
length.  Such  conditions  are  realized  in  a double-layer  system  when  the  splitting  energy 
Ar  is  small  in  comparison  with  the  Fermi  energy  eF.  Assuming  this,  we  derive  the 
boundary  condition  for  fPxPr(y): 


(ipy  + KP)fPxPy(yo)(-ipy  + KP+)  = ( ipy  + KP)fpx-py(yo)(iPy  + (5) 


where  pv  ~ hk+  ~ M_  is  the  absolute  value  of  the  electron  momentum  along  0y. 
For  a non-contacted  boundary  (edge),  the  diagonal  part  fPxPv(y)Pi  + fPxPv(y)Pr  of  the 
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distribution  function  obeys  the  boundary  conditions  fp'Pv(yo)  = fp'''-Pr(y o),  which  mean 
that  both  /-layer  and  r- layer  currents  through  the  edge  are  absent. 

3 Some  results  and  discussion 

The  boundary  condition  (5)  can  be  applied  to  various  problems  of  electron  transport 
in  double-layer  nanostructures.  To  demonstrate  the  effect  of  the  edges,  we  consider  a 
symmetric,  finite-width  ( —L/2  < y < L/2),  infinitely  long  (— oo  < x < oo  ) double- 
layer strip.  If  a weak  electric  field  E is  directed  along  the  strip,  the  linearized  kinetic 
equation  for  nonequilibrium  part  SfPxPv(y)  of  the  distribution  function  is 

1 [*»  s.Lp,(y)\  -~eE^s  & - (p?  + p2y)/2m) 

= - [( 1 + , $fp;Py O')]  + /2 T , (6) 

where  the  collision  integral  is  expressed  through  the  averaged  scattering  time  r and 
scattering  asymmetry  parameter  fi,  and  [...]±  denote  the  commutator  (-)  and  anti- 
commutator (+).  A solution  of  Eq.(6)  together  with  the  boundary  condition  (5)  allows 
to  find  the  resistivity  of  the  system.  The  edge  effects  dramatically  modify  the  resistivity 
of  the  strip  in  comparison  with  the  resistivity  of  a double-layer  plane  [5],  provided  that 
the  scattering  in  is  non-symmetric  (fj,  0),  and  qi  yt  qr.  The  most  important  qualita- 
tive modifications,  as  they  seen  from  Fig.  2,  are:  the  overall  suppression  and  asymmetry 
of  the  resistance  resonance  peak  [5],  and  the  oscillations  of  the  resistivity  upon  the 
background  of  the  peak.  These  phenomena  can  be  considered  as  manifestations  of  the 
quantum  size  effect,  which  occurs  in  double-layer  strips  at  L ~ Lp  = hvF/2T  (vF  is 
the  Fermi  velocity).  The  quantum  size  effect  is  suppressed  when  L exceeds  the  mean 
free  path  length  vft,  and  the  resistance  resonance  peak  becomes  the  same  as  for  the 
double-layer  plane. 

To  demonstrate  the  role  of  contacts,  we  consider  the  independently  contacted  double- 
layer systems  shown  in  Fig.  1 (the  layers  are  assumed  to  be  infinite  along  OX).  A solution 
of  the  kinetic  equation  together  with  the  boundary  conditions  (5)  with  complex  qi  and  qr 
allows  to  find  the  conductance  associated  with  the  double-layer  region.  The  conductance 
for  both  the  systems  has  oscillations  as  a function  of  A,  their  period  is  determined  by 
SAt  = 27 tHvf/L.  These  oscillations  are  of  the  same  origin  as  those  shown  in  Fig.  2, 
and  they  are  suppressed  at  L > vft. 


-L/2  y L!2  -L/2  LI  2 

Fig  1.  Two  kinds  of  double-layer  structures  realized  by  independent  contacting. 

In  conclusion,  we  have  derived  the  matrix  boundary  conditions  for  Wigner  distribu- 
tion function  of  electrons  in  tunnel-coupled  2D  layers,  and  demonstrated  the  quantum 
size  effects  occurring  in  double-layer  nanostructures  whose  dimensions  are  comparable 
with  tunnel  coherence  length  LT  = hvF/2T.  This  length  (Lr  ~ 100  nm  for  typical 
parameters  2T  ~1  meV  and  eF  ~ 10  meV)  is  large  in  comparison  with  the  quantum 
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Fig  2.  The  shape  of  the  resistance  resonance  for  several  values  of  2TL/Hvf  (the  values  are  shown 
near  proper  graphs)  calculated  in  conditions  (ITr/h)2  p1  and  vft  S>  L,  at  fi  = 0.8.  The 
dashed  line  shows  the  same  for  an  infinite  double-layer  plane.  The  inset  shows  the  model  of  the 
boundaries  used  to  estimate  qi  and  qr  (in  the  calculation  we  assumed  that  the  Fermi  momentum 
is  equal  to  4K/6L). 


wavelengths  of  the  electrons  on  the  Fermi  level.  The  effects  could  be  observed  ex- 
perimentally in  high-mobility  samples  at  low  temperatures,  when  the  electron  transport 
across  the  narrowest  dimension  is  close  to  ballistic. 
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